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Abstract 

We present an algorithm for the analytical evaluation of dimensionally regular- 
ized massless on-shell double box Feynman diagrams with arbitrary polynomials 
in numerators and general integer powers of propagators. Recurrence relations 
following from integration by parts are solved explicitly and any given double 
box diagram is expressed as a linear combination of two master double boxes 
and a family of simpler diagrams. The first master double box corresponds to 
all powers of the propagators equal to one and no numerators, and the second 
master double box differs from the first one by the second power of the middle 
propagator. By use of differential relations, the second master double box is 
expressed through the first one up to a similar linear combination of simpler 
double boxes so that the analytical evaluation of the first master double box 
provides explicit analytical results, in terms of polylogarithms Lia(— t/s), up 
to a = 4, and generalized polylogarithms Sa^bi~t/s), with a = 1,2 and 6 = 2, 
dependent on the Mandelstam variables s and t, for an arbitrary diagram under 
consideration. 
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1 Introduction 



The massless double box diagram shown in Fig. 1 with general numerators and integer 
powers of propagators is relevant to many important physical processes. The purpose 
of the present paper is to present an algorithm for the analytical evaluation of the 
general massless on-shell (i.e. for = 0, i = 1, 2, 3, 4) double box Feynman diagram 
in the framework of dimensional regularization [|I|, with the space-time dimension 
d = A-2e gularization parameter. 




k +pi 



l + Pi 



Pi 



l+Pi+Ps 



P3 



Figure 1: Planar double box diagram. 

The dimensionally regularized on-shell master double box Feynman integral (i.e. with 
all powers of propagators equal to one and no numerators) has been analytically 
evaluated in [§] in terms of polylogarithms Li^ (— t/s), up to a = 4, and generalized 
polylogarithms Sa,b{~t/s), with a = 1,2 and 6 = 2, dependent on the Mandelstam 
variables s = {pi + P2Y and t = (pi + Pa)^. In recurrence relations within 
the method of integration by parts (IBP) were explicitly solved and any given 
double box diagram was expressed as a linear combination of the master double box 
mentioned above, the second master double box which differs from the first one by the 
second power of the middle propagator, a family of double boxes with two contracted 
lines considered in shifted dimension, and vertex diagrams. 

In the present paper, we complete this algorithm by evaluating the second master 
double box and presenting crucial checks of our results by use of asymptotic expan- 
sions in the limits t/s —>■ and s/t 0. In the next section, we present definitions of 
the double box diagrams through integrals in loop momenta and a-parameters. Then 
we describe recurrence relations and their solutions in terms of the two master double 
boxes and a collection of simpler diagrams, in accordance with 0. In Section 3, we 
describe how the double boxes with two contracted lines are analytically evaluated in 
expansion in e up to a desired order. In Section 4, we present the analytical result of 
ref. for the first master double box and then use differential relations in order to 
express the second master double box through the first one up to a similar linear com- 
bination of simpler double boxes. This provides an explicit analytical result for the 
second master double box, in terms of the same class of functions as for the first one. 
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In Section 5, we use the general strategy of regions for expanding the double boxes in 
the limits t/s ^ and s/t — * 0. We present analytical algorithms for the evaluation 
of the hard-hard and collinear-collinear contributions to the expansion. We evaluate 
the LO and NLO (respectively order l/{sH) and l/(s^)) of the expansion of the first 
double box in the limit t/s ^ and find an agreement with our explicit result. In 
conclusion we discuss our results. 



2 Recurrence relations and their solution 

The general massless on-shell double box Feynman integral in ci-dimension can be 
written as 

K^''\P,ai,...,ar,s,t) 



(7r'^/2)2 (/2 + 2pi/)"i(/2 -2p20 
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where P is a polynomial, integers, k and / are respectively loop momenta of the 
left and the right box. Usual prescriptions, k"^ = k'^ + iO, —s = —s — iO, etc. are 
imphed. 

The a-representation of the double box with P = 1 is straightforwardly obtained 
(we omit in the following the polynom P and consider only scalar integrals): 



1 \ai+...+a7^ai + ...+a7+2-d 
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(2) 



where 



D = (ai + ^2 + a7)(a3 + "4 + tts) + tt6(ttl + "2 + "3 + "4 + "5 + tt?) ? (3) 
A = aia2(a3 + ^4 + ^5) + «3«4(ai + ^2 + a?) + «6(«i + «3)(«2 + ^4) • (4) 

To deal with Feynman integrals with numerators we use the fact that any polyno- 
mial in the numerators of the propagators can be represented as a differential operator 
with respect to some auxiliary parameters (see, e.g., [^]) acting on a scalar diagram. 
An outcome of this procedure is that any tensor integral is expressed in terms of scalar 
integrals but in different (shifted) space-time dimensions and with shifted indices of 
lines (see a detailed discussion in 0). This step is straightforward and formulae from 

can be easily programmed on computer. 

A more difficult part of the program is to express the so obtained double box 
integrals (in different dimensions and with all possible sets of indices) in terms of some 
master integrals and a family of simpler boundary integrals. It turns out that, in our 
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case, there arise only two master integrals K^f^ = K'''^\1, l,l,l,l,l,l;s,t), K^^2) = 
K^'^\l, 1, 1, 1, 1, 2, 1; s, t), and the boundary integrals are either vertex integrals that 
are evaluated in gamma functions and integrals with at least two reduced lines. 

Using the integration by parts method ^ the following relation can be obtained 
to reduce the index ai to unity: 

sail+ = 077+2" + a(i6^{2~ - 4") + ail+2" - {d - 2a2 - ai - a-j - oe) . (5) 

Hereafter we use the standard notation: j+ is the operator increasing/decreasing the 
index on the jth line by one unit, i.e. . . , a^, . . .) = K{. . . , ± 1, . . .). 

Three similar relations obtained by permutations of lines can be used to reduce 
indices of lines 1,2,3 and 4 to one. Next we can reduce indices of lines 5 and 7 with 
the help of the following relations: 

(d — 2 — 2a5 — a4 — 03)055'*' = {d — 2 — 2aQ — — 03)066'*' 

+ (05 — 05)044+ + (05 — 05)033+ + 04054+6+2^ + osOgS+G+l" , (6) 

[d — 2 — 207 — 02 — 01)077+ = [d — 2 — 2o6 — 02 — oi)o66+ 

+ (07 — 05)022+ + (07 — 06)oil+ + 02052+6+4^ + oiOgl+G+B". (7) 

Using the above recurrence relations we can bring indices of lines 1,2,3,4,5,7 all to 
unity so that only Og can be greater than one. 

Our relation to reduce the index of hne 6 reads 

t{d-Q- 2o6)(o6 + 1)056++ = 

— {d — 5 — a(i)[3d — 14 — 2o6 + 2o6-]o66+ 
2 

+ -(d — 4 — oe) (d — 5 — oe) 

s 

-(d-A- a^Md - 5 - oe) + 2-006+ 

s s 

2t(o6 + l)o66++ + 2(c/ - 4 - 06)o66+]3+|l- 

+ (d-6)7-d-, (8) 

where d~ decreases the dimension of space-time by 2. Note that the dimension can 
be effectively shifted only by an even integer number (see detailed discussion in 0). 
We stress that this formula is valid only if the indices of lines 1,2,3,4,5,7 were already 
reduced to unity. Note that in the left-hand side of (^) there is 6++, rather than 6+. 
This means that the index of line 6 cannot be always reduced to one but generally 
to one or two. One can also get rid of d~ in the above formulae replacing it by 
(1+ + 2+ + 6+)(3+ + 4+ + 5+) + 6+(l+ + 2+). 
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Let us comment shortly on how and other similar relations can be derived. 
One can start from the integral with numerator 2kp2- Since 2kp2 = {k — ^2)^ ~ k"^ we 
can eliminate this numerator by canceling lines 4 or 5 (see Fig. 1), and the resulting 
integrals are simple. On the other hand one can use the machinery mentioned after 
(^) (see 1^, ^) which expresses tensor integrals as differential operator acting on scalar 
integral in some dimension d + 2n. It is more convinient in our case to start from the 
dimension d — 2. Then we have 



(9) 



where dj = d/drn^ takes the derivative with respect to the square of the mass on 
jth line. (After differentiaton all masses are put to zero). In the right-hand side of 
the above formula, there are scalar integrals with increased indices. Therefore we can 
apply reduction formulae (|§)-(0). The resulting relation will involve terms like 6"'""'" 
and, after some transformation, one can come to (P). 

To complete the reduction procedure we should bring the two master integrals 
(which can appear in shifted dimensions) to the integrals in the generic d = 4 — 2e 
dimension. Thus we need also a relation that reduces the dimension of the space-time. 
This can be obtained by inverting the identity 



ai. 



D{d)K^'^\a^, 



(10) 



where K^'^'^ is defined in (|Tp, D{a) is given by (|^), and d denotes a family of dif- 
ferential operators acting on auxiliary masses of the lines: di = d/dmf. (After the 
differentiation, all these masses are put to zero.) 

(p!0[) is valid for arbitrary indicies Oj and can be derived from the a-representation 
of the Feynman integral (see e.g. This relation however increases the dimension 
by 2 units. To find a relation decreasing dimension we have to compute operator 
which is inverse to D. These can be done with the help of already listed above 
reduction formulae. Instead of giving general form of inverse to (|T0|) it is enough to 
give it for the master integrals. Since we have two master integrals there are two 
relations ISfl 



^(i)-A 
^(2) - A 



-021 (K, 



^(1) 

(d-2) 
(1) 



where operators /j are given by 

A'' = 

2+4+ - 




I T^(d-2) 



^(1); - "i2v-'>'(2) 



Ad) T^{d)y 
J 2 -"-(1); 



(d-2) 
(2) 



J 2 



2+6+ + 4+6+ + 4+7+ + 3+7^ 



(11) 
(12) 



4 



+-(2+5+ + 5+6+ + 5+7+) - -^(d - 5)(3s + 2^)(2^ 
s sH 

2 



3+6+7+ - g^^/ g) (^^(^ ~ ^) + ^^^^ ~ 14))3+6+|l- 
+ ^7-d-, (13) 



d-6 
3 



/2 



+ (-(2+3+ + 2+4+ + 3+7+ + 4+7+)6+ + -(2+ + 4+)6^ 

ys s 

+ ^^^j~^f (2+ + 7+ + 2-6+)5+6+ 



2(d-5)(d- 7) 



(s(3rf - 20) + 2t{d - 6)) (2+ + 7+)6^ 



sH{d-Q){d-S 
^ ^ ^ + 3^)7^6^^ + + - 13). 



+ {|^«' - + «(2.- 13)))7-d- . (14) 

aii = J^(c/-5)2(3s + 2t), (15) 

ai2 = --(4rf-21) - -(3c/- 16), (16) 
s t 



^21 



(d-5)2(d-7)/8(2d-13) 6(3rf - 20) 



^ o^za - loj ^ o^oa-zuj s^ ^^^^ 



st{d-8) ^ s t 



1 T 

022 = -^TTTi 7 (Ss^iSd - 16) (3d - 20) + 6st(5d^ - 59d + 172) 

sH^[d — 8) ^ 

+4t2(d-5)(rf-6)), (18) 
16(s + t)(c/-5)^(c/-6)(ci-7) 

;W^8^ • ^ ^ 

Formulae (^-(|19D solve the problem of the reduction of any planar double box to 

two master integrals Kf^i-^ , K(2) and a set of simpler integrals with reduced lines. From 

now on we omit the superscript {d) because we shall deal with integrals in d = 4 — 2e 

(non-shifted) dimensions. 

We have checked the reduction scheme (p|)-(|T9|) by expanding the integrand in 

two regimes, s/t — >• and t/s — > 0, and evaluating the resulting one-scale integrals 

by the method described in Sect. 5. 
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3 Double boxes with two contracted lines 



In the framework of the reduction procedure presented in the previous section, the 
boundary values for general double boxes are either two master double boxes, or 
vertex diagrams, i.e. at as = or aj = 0, or double boxes with two contracted lines. 
The latter can be of the following two types: ai = = (or the symmetrical variant 
0^2 = os = 0) shown in Fig. 2, or = = (or the symmetrical variant ai = a2 = 0) 
shown in Fig. 3. Let us call them respectively the box with a diagonal and the box 
with a one-loop insertion. Note that they generally arise in a shifted dimension. 

We should mention that there are no integrals with only one contracted line left. 
As far as one of the lines 1,2,3 or 4 is contracted one can still proceed with reduction 
formulae from Sect. 2 or apply the standard "rule of triangle" to reduce one more 
line. Thus the two master double boxes plus boxes with two contracted line form the 
basis of integrals. 
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Figure 2: A box with a diagonal. 




Figure 3: A box with a one-loop insertion. 



These cases without two lines are much simpler than the master double boxes. 
Using a-representation (0) and representing one of the functions involved into the 
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Mellin-Barnes (MB) integral 

1 1 1 /•+'^ , Y 



— / dw——r(u + w)r(-w), (20) 



{X + YY r(i/) 

we obtain the following results: 

K^'^\0,a2,a3,0,a5,ae,a7]d+n;s,t) = ]-[ p(Q.^(_g-)a-n-4+2e 

r(2 - ag - as - e + n/2)r(2 - - e + n/2)T{2 - 02 - ay - e + n/2) 



X: 



r(6 — a + 3n/2 — 3e)r(4 — — — ae + n — 2e)r(4 — a2 — ae — ay + n — 2e) 
X — / dz (i/s)T(a5 + z)T{a7 + z)T{a - 4 - n + 2e + -z)r(-z) 

2Tn J-ioo 

xr(4 — a2 — as — ae — ay + n — 2e — -2)r(4 — az — a^, — — a-j -\- n — 2e — z) (21) 
and 

i2(-l)" 

ii:(ai,a2,0,0,a5,a6,a7;rf + n;s,t) = Y{V{a){-sY~'^~'^'^'^'- 
r(2 - ag - e + n/2)r(2 - ag - e + n/2) 



X 

27ri 



r(6 - a + 3n/2 - 3e)r(4 - ag - ae + n - 2e) 

\ r+ioo 

— dz {t/syv{a7 + z)T{-A + a - n + 2e + ^) 

TTZ J —ioo 

xr(— 2 + as + ae + e — n/2 + z)r{A — ai — — qq — ar + n — 2€ — z) 

xr(4 - a2 - as - ae - ay + n - 2e - z)r(-z) , (22) 

where a = J2iCi-i- The contour of integration is chosen in the standard way: the poles 
with the r(. . . + 2;)-dependence are to the left of the contour and the poles with the 
r(. . . — 2;)-dependence are to the right of it. 

Then each of these boundary integrals with given values of integer indices is de- 
composed into 'singular' and 'regular' parts. For the diagonal crossed boxes, the 
singular part is written as minus the sum of the residua of the integrand at the points 
j — 2e, with j = — max{a2, as} — as — ae — ay + 4 + n, . . . , — 1, plus the sum of the 
residua of the integrand at the points j — 2e for j = 0, . . . , 4 + n — a2 — as — as — ag — ay . 
For the boxes with one-loop insertions, the singular part is written as minus the sum 
of the residua of the integrand at the points j — 2e, with j = — max{ai, 02} — 05 — 
ae — ay + 4 + n, . . . , — 1, plus the sum of the residua of the integrand at the points 
J — 2e for j = 0, . . . , 4 -|- n — ai — a2 — as — ae — ay. 

The regular parts are given by MB integrals where no gluing of poles of gamma 
functions with +z and —z dependence arises. They can be written as MB integrals for 
— 1 < Kc{z) < with an integrand expanded in a Laurent series in e up to a desired 
order. Then these integrals are straightforwardly evaluated by closing the contour of 
integration to the right and taking residua at the points z = 0,1,2, .. .. At this step, 
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one can use a collection of formulae for summing up series presented in ref. . The 
evaluation of both the singular and the regular parts is easily realized on computer. 

Note that in the recurrence relations of the previous sections, these boundary 
double boxes can arise with coefficients that involve poles up to the second order in e 
so that the expansion up to is here necessary. However the 'master' boxes with a 
diagonal or one-loop insertion enter with coefficients finite in e so that it is sufficient 
to compute then, in expansion in e, up to the finite part. Let us, for example, present 
an analytical result for the master box with a diagonal in d dimensions: 

i2 



Ue 



i^(0, 1, 1, 0, 1, 1, 1; d- s, t) = ^—-^Kois, t, e) , (23) 



where 



+ 



2e2 

2Li3 (-t/s) - 2 ln(t/s)Li2 (-t/s) - (\n\t/s) + vr^) ln(l + t/s) 

1 



+- ln%t/s) + In(-s) \n\t/s) + In(-t) - 2C(3) 
3 

+4 {S2,2{-t/s) - Ht/s)Si,2{-t/s)) - 4Li4 i-t/s) 
+4 (ln(l + t/s) - H-s)) Lis i-t/s) 

+2 [Hit/s) + 2 In(-s) ln(t/s) - 2 ln(t/s) ln(l + t/s)) U2 (-t/s) 
+2 Q ln='(t/s) + In(-s) \n\t/s) + vr^ In(-t) - 2C(3)) ln(l + t/s) 

- (ln2(t/s) + TT^) ln2(l + t/s) - - ln^(t/s) - - In(-s) ln='(t/s) 

2 3 

\n\-s) + —nA Hit/s) - TT^ \n^i-s) - 2n^ In(-s) ln(t/s) 
X2 / 

+4C(3) In(-t) - ^ . (24) 

Here Li^ (z) is the polylogarithm and 

S (z) - r ^^^-\t)\n\l-zt) 
^^^'^'^ - ia-l)\b\ Jo 't 

a generalized polylogarithm Using known formulae that relate polylogarithms 
and generalized polylogarithms with arguments z and 1/z [j^, one can rewrite this 
and similar results for the master double boxes in terms of the same class of functions 
depending on the inverse ratio s/t. 

We do not give explicit result for the box with a 1-loop insertion since it is certainly 
a simpler object. Indeed the 1-loop onsertion can be trivially integrated and we have 
a 1-loop box where one of the indicies is equal to e. Such kind of integrals can be 
expressed using standard methods in terms of the hypergeometric function 2F1 (see 
e.g. 
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4 Master double boxes 



The first master double box 



Ue 



-7e<; 



,2 



4 



K{1, . . . , 1; s, t) = ^_^^2+2e(_^) -^i(Vg> e) , (26) 

has been evaluated in ref. by use of a-parameters and resolving singularities in a 
5-fold MB integral: 

/ N 4 51nx / , 9 5 n\ 1 

i^i(x,e) = -- + -^-(^21n2x--vr^j- 

3 11 21 65^,^, \ 1 4,4 ^ 2, 2 88^,^,, 29 

- - In^ X H vr^ In x ((3) - + - In^ x + Gyr^ In^ x ((3) Inx H vr 

V3 2 3 7 e 3 3 ^ 30 

2Li3 (-x) - 2 InxLis (-x) - (in^ x + vr^) ln(l + x) 

-4(5'2,2(-a;) - lnx5'i,2(-a;)) + 44Li4 (-x) - 4 (ln(l + x) + 6 Inx) Lis (-x) 

+2 {\-n^ X + 2 Inx ln(l + x) + y ^r^) Li2 (-x) 

2 

+ (in^ X + vr^) ln2(l + x) - - (4 In^ x + Svr^ In x - 6C(3)) ln(l + x) . (27) 
To evaluate the second master double box, i.e. 

K{1, 1, 1, 1, 1, 2, 1; d- s, t) = -^—^^K^it/s, e) , (28) 

let us take first derivatives in t of the two master double boxes. Using a-representation 
(|^) we obtain 

^ir(l,...,l;ci;s,t) = -ir(l,l,l,l,2,2,2;ci+2;s,t), (29) 

^ir(l,l,l,l,l,2,l;rf;s,t) = -2if(l,l,l,l,2,3,2;rf + 2;s,t). (30) 

We now use the results presented in Sect. 2 to express both right-hand sides as linear 
combinations of the two master double boxes, vertex diagrams (two in the first case 
and three in the second case) and a numerous family (around fifty terms in each case) 
of diagonal crossed boxes and boxes with one-loop insertions. Substituting explicit 
result (^) into the first equation and evaluating all the terms in the right hand side 
as explained in Section 3 we obtain an analytical result for the second master double 
box: 

K2{x,e) = ^-5(lnx-2)^+ (^21n2x-141nx-^(7r2 + 4)) ^ 
+ (^^ln^x + 81n2x+ (^yvr^ + 14^ inx - 2 - 37r2 - ^((3)^ ^ 



4 / 88 \ 

-- ln=^ a;(ln x + 1) - 2 (Stt^ + 4) In^ ^ + M + Qtt^ + y C(3) j In x 



+20 + 127r2 - — vr^ + -((3) 



+x 
1 

+ 6 

+ 

-4x 



7 1 
— r + (81na; - 33) — + (261na; + 6 + 



21 



-TT 



-32 In^ X - 4(21 + 267?^) In x + 180 + 2mTx'^ + 904C(3) 



2Li3 (-x) - 2 lnxLi2 (-x) - (in^ x + tt^) ln(l + a 
8 (Lig (-x) - InxLis (-x)) - 4 (in^ x + vr^) ln(l + 



X 



+4 {S2,2{-x) - lnxSi,2(-x)) - 44Li4 (-x) + 4 (ln(l + x) + 6 Inx - 2) Lig 



-X 



2 ( In^x + 21nxln(l + x) -41nx+ — tt^ ) Lig 
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-X 



In^x + TT^ ln^(l + x) 



+ ( ^ In^ X + 4 In^ X + y vr^ In x + An'^ - 4C(3) ) ln(l + x) . 



(31) 



Proceeding in the same way with the second recurrence relation (|30D, and inserting 
there our analytical results for the master double boxes we eventually obtain an 
identity of the left-hand and the right-hand sides. This fact turns out to be a very non- 
trivial check of the recurrence relations, their solutions and our analytical expressions 
for the master double boxes. 



5 Asymptotic expansions of the double box 

We still want other checks and are going to compare our results with what can be 
obtained by expanding the first master double box in various limits. To expand the 
double box diagrams in the limit t ^ let us use the strategy of regions: 

(^) Consider various regions of the loop momenta and expand, in every region, the 
integrand in a Taylor series with respect to the parameters that are considered small 
in the given region; 

(m) Integrate the integrand expanded, in every region in its own way, over the 
whole integration domain in the loop momenta; 
{Hi) Put to zero any scaleless integral. 

In the off-shell and off-threshold limits, this strategy leads to the well-known 
explicit prescriptions fTl]] (see a brief review [|12|) based on the strategy of subgraphs. 



Although the strategy of subgraphs was successfully applied to some on-shell limits 



13, H, the strategy of regions looks generally more flexible. In particular, it proved 



to be adequate for constructing the threshold expansion [|l^ . 



Let us choose, for convenience, the external momenta as follows: 



Pi,2 = (+g/2, Q/2, 0, 0), r = + p3 = (T/Q, 0, JT + T^/Q\ 0) , (32) 
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where s = —Q^ and t = —T. The given hmit is closely related to the Sudakov limit 
so that the following standard regions happen to be typical for it: 



hard (h) 

1- collinear (Ic) 

2- collinear (2c) 
ultrasoft (us) 



k+r^Q, k^r^T/Q, kr^VT, 

k. r^Q, k+^ T/Q, k^Vf, 
k^T . 



Here k± = ko ± ki, k={k2, k^). We mean by A; ~ Q, etc. that any component of k^ 
is of order Q. 

It turns out that the (h-h), (Ic-lc) and (2c-2c) are the only non-zero contribu- 
tions to the asymptotic expansion in the limit t/s —>■ 0. In particular, all the (c-h) 
contributions and all the contributions with ultrasoft momenta are zero because they 
generate scaleless integrals. 

The (h-h) region generates the contribution given by Taylor expansion of the 
integrand in the vector r. Every diagram from this contribution corresponds to the 
forward scattering configuration, ps = —pi and p^ = —p2, and can be evaluated 
for general e in gamma functions by resolving recurrence relations following from 
integration by parts The first step of this procedure is to reduce an index 05 
or 07 to zero and thereby obtain vertex massless diagrams. The latter reduction, in 



the scalar case, was constructed in ref. fl^. (In the case without numerators, the 



reduction of the forward scattering double boxes was presented in ref. p!6[] .) 

We have constructed two different procedures for the evaluation of the (h-h) part: 
along the lines of this standard recursion and also by expanding the integrand of 
the a-representation in the variable t and using tricks with shifting dimension. We 
have implemented both methods and checked that they give the same results for first 
several coefficients. 

We describe now how one can expand the integrals in the variable t. The most 



suitable method to achieve this is the one proposed in |T^. For the expansion in t/s, 
we have 



00 



irW(s,t)^='^Ei- (--)'Qiir('^+2^)(s,0), (33) 

where Qt is the differential operator acting on the masses of the lines. If we denote 
d/dmf as di then 

Qt = d^d^dj. (34) 

After the differentiation in (|53|), all these auxiliary masses are put to zero. As a result 
each coefficient in the expansion ( p3D consists of integrals depending only on one scale 
s. These can already be evaluated analytically by the standard "triangle" rule, i.e. 
using the integration by parts 0]. 
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The integrals K^^~^'^^\s, 0) belong to the class of primitives, i.e. they can be evalu- 
ated in terms of F-functions. The problem here is that the repeated application of the 
triangle rule brings more and more powers of 1/e and therefore a deeper expansion of 
F-functions is required. To keep things under control one can use e.g. an algorithm 
described below. With the help of this algorithm, the depth of the e-expansion is 
kept at the level of 6. The reduction proceeds as follows: 

(i) Use the relation 

{2d — 2a^ — 2aQ — 2a'j — ai — 02 — eta — 04) 
= + a22+7" + agS+S" + 044+5" . (35) 

(This is the same first step as in the above mentioned standard recursive procedure.) 
With its help, we can get rid of either line 5 or 7 and thereby reduce the double box 
to a planar vertex. Relation (|35|) has the feature that the left-hand side passes only 
once through the "critical point" (when the expression in parentheses is proportional 
to e). Therefore at most a single power oil/e will be generated in the course of this 
step of the recursion. 

(ii) Reduce indices of the lines 3 and 4 to one with the help of 

5038+ = 055+4- + a66+(4- - 2") + 038+4- -{d- 204 - 03 - 05 - oe) , (36) 
5044+ = 055+3- + a66+(3^ — 1^) + 044+3- — {d — 203 — 04 — 05 — og) , (37) 

This step brings no new powers of 1/e. 

(iii) Shrink the line 4 (or 3) using the triangle rule 

{d - 2aQ - oi - 02) = ail+(6- - 3") + a22+(6" - 4") . (38) 

Here at most one additional power of 1/e is generated. 

(iv) The index of the line 1 is then reduced to unity by 

ail+ = (rf - 1 - 2a7 - ag - fli - 02)^+ + 056+7+5- - 032+ - 056+ . (39) 

No new powers of 1/e arise here. Note that the line 7 could be absent before this 
step but now it appears again due to the 7+ term. 

(v) Apply the triangle rule 

{d - 2a3 - as - ag) = a66+(3- - 1") + 055+3- (40) 

to shrink either line 1 or 3. Note that a dangerous term 6+1- which leads to the 
oscillation of the left-hand side around e is harmless here because the index of the 
line 1 is already equal to unity due to (iv). 

(vi) Now we have to evaluate the diagram without lines 1 and 4. (Other cases are 
trivial). To simplify it we use (^) which now looks like 

{2d - 2a5 - 2aQ - 2a7 - 03 - 03) = a22+7- + 038+5- . (41) 
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This completes the algorithm. 

Thus steps (i),(iii) and (vi) in the worst case bring a factor 1/e each. Other three 
powers of 1/e arise when evaluating primitive integrals. Therefore the total depth of 
the e-expansion of the F-functions must be not greater than six. Note that diagrams 
without line 6 must be separated from the very beginning and evaluated immediately 
without any recursion in terms of F-functions. The point is that when evaluating 
primitive diagrams with og = we obtain poles of the fourth order instead of usual 
third order poles. 

Let us now describe how an arbitrary term of the (c-c) contribution to the ex- 
pansion of the (first) master double box can be evaluated. If we consider k and / 
1-collinear then fc^, and P2I are of order T while pik and pil are of order Q^. 

Moreover, (/ + r)^ = + 2/r — T ~ (Z + f)^, where r = pi + p^ and 

f =(T/(2Q),-T/(2Q),Vt,0), (42) 

with 2pif = 0, 2p2r = = — T. Thus the (Ic-lc) contribution is obtained by 
expanding propagators 1/(A;^ + 2pik) and l/(/^ + 2pil) in Taylor series respectively in 
k'^ and Z^, and by expansion also in Taylor series in 2pir. (Note that we are dealing 
with a function of three kinematical variables, 2pir, 2p2r and r^. So, we expand the 
integrand (e.g. in the a-representation) in 2pir and then put 2pir = T.) Actually 
we want to expand the master box with all Oj equal to 1. Therefore only the leading 
term in the Taylor expansion in k"^ is non-zero because, starting from the next order, 
the factor A;^ cancels the propagator l/k"^ and we obtain a zero scaleless integral. 

In the following it is more convenient to consider integrals with general arbitrary 
indicies Oj. Thus, we need integrals of the type 

, , 1 f d V r r d'^kd'^l 

J{ai,...,aQ;d,s,t) = — r tttp 



X 



asiydXj J J (7r'^/2)2(-2pi/)«i(-/2 + 2P2O''' 

(/2)«9 



(-2piA;)«3(-A;2 + 2p2kY^{-k^)^^'{-{k - /)2)"6(_(/ _ r)2)ar 



, (43) 



x=o 



with X = 2pir, for integer a^. However this integral taken alone is not generally 
regularized dimensionally. Only if we add the corresponding symmetrical contribution 
(i.e. for oi ^ 02, as 04) we shall have a result that exists within dimensional 
regularization. An efficient way to deal with this problem is to introduce an auxiliary 
analytical regularization which allows to consider the above terms separately. Let 
us introduce it into the lines 1 and 2 (although we could choose as well 3 and 4, 
or all four of these lines), i.e. with ai ^ ai + Xi,a2 — ^ 02 + X2 plus a symmetrical 
contribution which is given by interchanging ai + Xi and 02 + X2. Only in the sum 
we may switch off this regularization, i.e. let Xi — >■ 0,X2 — > 0. 
The general integral (^3]) at ag = ag = can be represented as 



^ ^ , N i\-l)''T(a' -d) 

J{ai, . . . ,aj,{J,U;d,s,t) — 



n.^l,3r(a,)(Q2)ai+a3Ta'- 
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X /da'5 E«^-l n «r"'(«5«6«7)'-'^'5'^-^"^r^/^3"% (44) 



where 



D = (^2 + a7)(a4 + as) + «6(«2 + "4 + "5 + «?) , (45) 
Di = a2(a4 + tts) + a6(a2 + "4) , (46) 
D3 = ai{a2 + a?) + a6(«2 + "4) • (47) 

(48) 

Moreover, the integral above is represented in terms of parameters a' = , i = 
2,4, 5, 6, 7}, a = Y^^i, 

a' = J2ij^i,3 di ■ In the argument of the delta function one can put the sum of an 
arbitrary subset of ai ,i = 2, 4, 5, 6, 7. 

The integrals with ag > are obtained by the replacement 

where Bi = a'T^a^ + ^5 + as) and -B3 = a^aj, subsequent expansion in Taylor series 
in X and keeping terms of order 

Starting from ( PD and it generalization for ag > we can arrive at a triple MB 
integral in the following way: 

1. Choose the delta function as 5(^4 + aQ — l). 2. Represent 1/(D3)"='"'"'" in a MB 
representation, using the decomposition D3 = Dq + a^a-j. 3. Integrate in a-j. Note 
that the function Di + a^a^ arises. 4. Represent {Di + a^a^)'" via MB representation. 
5. Represent {Di = Dq + 020:5)'" via MB representation. 6. Integrate in and then 
in a2- 7. Integrate in ^ = 04 with = 1 — C,- 

Then we obtain a triple MB integral of a ratio of F-functions. This integral 
is evaluated, in expansion in e, by the standard technique of shifting contours and 
expansion in MB integrals. First, singularities in xi — X2 are localized. Second, the 
same technique is applied for picking up singularities in e. As a result we end up with 
a collection of explicit terms plus integrals which are finite in e. These last integrals 
can be expanded in e up to the desired order. So, in the end, various integrals with 
F, -0 and their derivatives are evaluated (see examples in 0). 

The integrals with ag > can be represented in a-parameters in a cumbersome 
way. Using IBP 0] (starting from the integral of {d/dk) k) we obtain the following 
recurrence relation: 

(c/ - as - 04 - 2a5 - - a44+5" - a66+(5" - 9+)) J = . (49) 

Note that our general integral is zero when at least one of parameters (I5, Og, CLj IS a 
non-positive integer so that this relation is of no use in its primary form. However we 
may turn to a new integral which is defined without F(a6) in the denominator, i.e. put 
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J(. . . a6, . . .) = J'i.- ■ ■ ctS, . . .)/T{aQ). Then ( ^9[ ) applied to the new integrals J' can be 
used to decrease the number ag to zero. We shall generally deal with non-positive 
after that. 

To evaluate the (c-c) contribution to the expansion of the first master double 
box it suffices to consider CL2 — CLs — (34 — O5 — fly — 1, and general integer 
ai,ae and as,aQ > 0. Then ag can be reduced to zero using (|49|). In particular, 
in the leading order of the expansion we meet J(l, . . . , 1, 0, 0) and, in the NLO, 
J(2, 1 . . . , 1, 0, 1) and J(l, 1 . . . , 1, 1, 0). Using (|9|) we then express J(2, 1 . . . , 1, 0, 1) 
through J(2, 1, 1, 1, 1, 0, 1, 0, 0) (defined at negative ag as explained above). 

We have evaluated the LO and NLO of the (c-c) contribution (and added the LO 
order of the (h-h) contribution which is really NLO for the whole expansion) of the 
master double box integral in the limit t/s —>■ 0. After that we have found complete 



agreement with the first two terms of the expansion of the explicit result (27): 



4 51nx /,o 5n\l 



3 11 21 65^,^, \ 1 
-(-ln3. + -.Mn.--C(3))- 

4 88 29 

+ - In^ X + Qtt^ In^ x C(3) Ina; + — vr^ 

3 3 ^ ^ ^ 30 

+2x - fln^x - 21nx vr^ + 2 

.e ^ 

(4 In^ X + 3 In^ X + (Svr^ - 36) Inx + 2(33 + Svr^ - 3C(3))) 
+0{xHn^x). (50) 

The expansion in the limit s/t — *• has the same structure as in the previous case: 
only (h-h) and (c-c) contributions are non-zero. However, in this case, there are three 
(c-c) contributions, and the poles with respect to an auxiliary parameter of analytic 
regularization happen to be up to the second order (as in the case of the non-planar 
vertex diagram in the Sudakov limit — see ||ir 



6 Conclusion 

The on-shell double box provides a curious example where the analytical evaluation 
is simpler than the evaluation by means of asymptotic expansions in some limits 
(in this case, t/s and s/t — > 0). We have met rather inconvenient recurrence 
relations for diagrams that enter the coUinear-coUinear contribution. Even the global 
recurrence relations for the unexpanded diagram turned out to be simpler, with the 
solutions described in Section 2. Although a systematical evaluation of the (c-c) 
contribution ot an arbitrary order in the expansion is certainly possible, it would 
be difficult to guess the analytic form ot the result by evaluating and studying first 
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terms of the expansion. Still we really used, in two points, the method of asymptotic 
expansions for crucial checks of our explicit procedure. Firstly, we have used the 
possibility to evaluate an arbitrary term of the (h-h) contribution (both in the limits 
(t/s — i> and s/t ^ 0) for checking (global) recurrence relations and their solutions 
presented in Section 2. This was possible because any recurrence relation derived 
within integration by parts commutes with the (h-h) expansion which, by definition, 
is performed under the sign of the integrals involved (either in the integrals in the 
loop momenta or in the a-parameters) . Secondly, we have evaluated the LO and NLO 
contributions of the first master double box and successfully compared them with the 
expansion of the analytical result. 

In fact, we could avoid the rather cumbersome evaluation of the first master 
double box by starting from equations ( pOf ) and (0), applying solutions of recurrence 
relations of Section 2, evaluating all simpler diagrams, eliminating the second master 
double box by use of the second equation and arriving at a second order differential 
equation for the first master double box. It would be possible to solve this equation 
by expanding in t/s. But then we would need boundary conditions to solve it. Here 
we could insert the first two orders of the expansion in the limit t/s evaluated by 
means of the strategy of regions as explained in Section 5 (where we have met only 
3-fold rather 5-fold MB integrals.) However, in this case, we would not have crucial 
checks for the obtained results. 
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